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An analytical study of phase contrast imaging is performed. It is 
shown that the complex disturbance in the image plane can be repre- 
sented as a convolution of the object disturbance and the Fourier 
transform of the transmission function of the phase plate. The simple 
theory of the phase contrast microscope is then derived as a limiting 
case of this more general result that is applicable when the size of the 
phase object is small compared to the area of the entrance aperture of 
the system. The response of a general system to several simple large 
phase objects is also examined, and it is shown that qualitative infor- 
mation about these objects can be obtained from the intensity pattern 
when the phase perturbations are small, providing the background is 
sufficiently uniform, and the size of the phase spot on the phase plate is 
carefully chosen. The study provides insight into the type of perform- 
ance that can be achieved, for example, if phase contrast imaging is 
used to extract phase information from an optical memory or if it 
is used as an experimental tool to study the qualitative behavior of 
such phenomena as clear air turbulence. 

I. INTRODUCTION 

One of the most popular techniques for converting a spatial phase 
variation into a spatial intensity variation is the phase contrast 
imaging scheme originally proposed by F. Zernike 1-3 in 1935. Pro- 
ceeding mainly from a heuristic point of view, Zernike recognized 
that, with a conventional imaging system, most of the direct light 
passes through a small region R in the focal plane while, in many 
instances, most of the diffracted light is scattered away from this 
region. If a plate whose transmission function is given by 

T = \ae ia in R .^ 

i 1 otherwise 
1943 
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is placed in the focal plane, the direct light can be modified and made 
to interfere with the diffracted light to produce an intensity pattern 
in the image plane. Moreover, if the phase perturbations $ that are 
introduced by the object are small enough so that the approximation 

e'* « 1 + to 
holds, the intensity in the image plane and the phase of the object 
may be linearly related. These concepts led to the design of the phase 
contrast microscope. 

The key assumption in the theoiy is that the effect of the phase 
plate on the diffracted light is negligible. In 1953, H. H. Hopkins 4 
rigorously demonstrated that, if this assumption holds, a phase object 
will produce an intensity pattern given by 

I = ± 2 \ae ia - 1+«'*| 8 . (2) 

In this expression, M is the magnification of the system, and it is 
assumed in the derivation that the object is illuminated by a unit 
amplitude monochromatic plane wave. When </> is small, a linear 
relationship, 

1 = ]|p(a 2 + 2^sina), (3) 

results, as predicted by Zernike. f 

There are many potential applications of phase contrast imaging 
in which the basic assumption employed in Refs. 1-5 will not be 
valid. Therefore, it is interesting to examine what results will follow 
if these assumptions are not made. In this regard, Hopkins rigorously 
studied the diffraction images of circular disks under coherent illumi- 
nation, and M. De and S. C. Som, 7 and De and P. K. Mondal, 8 in- 
vestigated the images produced by similar objects under incoherent 
illumination. In this paper, we will derive the theory from the point of 
view of Fourier optics and demonstrate that, in general, the complex 
disturbance over the image plane can be represented by a convolution 
of the object disturbance and the Fourier transform of the transmission 
function of the phase plate. The approximate theory, equations (2) 
and (3), is then shown to be a limiting case of the general theory, 
applicable if the size of the phase object is small compared to the 
entrance aperture of the system. Such is usually the case when the 



t A related analysis is also given in Born and Wolf, 5 Chapter 8. 
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magnification of the system is large. For larger objects, the convolution 
integral yields a complicated relationship between the phase of the 
object and the intensity of the image. 

The response of the system to several simple large phase objects is 
then studied, and it is shown that qualitative information about the 
phase distributions of these objects can be obtained when the phase 
perturbations are small if an appropriate phase plate is employed. 

The work that is discussed in this paper complements the experi- 
mental investigations carried on at Bell Telephone Laboratories by 
N. J. Kolettis. 

II. GENERAL THEORY 

The optical system that we will consider is shown in Fig. 1. A 
circular entrance aperture with radius R is located a distance d in 
front of the lens, and a phase plate with a transmission function 



T f (r f ) = 



ae 



r, ^R" f 



(4) 



R°r <r f ^ R !o 
10 R l0 < r, 

is centered in the focal plane. The phase object is assumed to occupy 
some portion of the entrance aperture, and the illuminating light is 
taken to be a normally incident unit amplitude monochromatic plane 
wave. 

Let the complex disturbance in the object plane be denoted by 
U (x , y ) where U„ is assumed to be zero at points outside the 
entrance aperture. Then, from the Fresnel approximation, it follows 



<o ,__d - + f * » d 
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Fig. 1 — System geometry. 
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that the complex disturbance in the focal plane to the left of the 
phase plate is 

N exp [jk(d + / + n A )] [jk (. d\ 2 ,."] 

U,(x, , y f ) = H ix/ ; — exp |jj [1 - j)(x f + y r ) J 

• jj f/„(x„ , 2/ )Pl(» o + j x, , y + j i/,J 

• exp -^ (a^r, + y y f ) J dz d?/ . (5) 

In this expression / is the focal length of the lens, A is the wavelength 
of the illuminating light, k is the wavenumber 2tt/A, and /cnA is the 
phase shift introduced by the lens. P L is the pupil function of the lens 
which assumes the value of unity over the region covered by the lens, 
and is zero elsewhere. 

Now, the complex disturbance on the right side of the phase plate 
is simply U f T f , and, therefore, if the Fresnel approximation is also 
used to describe the propagation from this plane to the image plane, 
the complex disturbance in the image plane becomes 

Ufa , yd 

= Wd /£ /£ **■ ■ ^ + 1 »" * + 1 y ') Tfix ' ' yr) 

•exp |-^ [(Mx + xdx, + (Mj/ + y t )y,]\ dx dy dx, dy, , (6) 
where 

K = -exp [ifc(d„ + / + d + n A )] exp |j| (.x- + t/?)J- 
Here we are assuming that the distances d , d, and / obey the lens law, 

and we have denoted the magnification of the system, d/f, by M. 

The integral in (6) is complicated by the coupling of the coordinates 
x , x f and y , y f via the pupil function of the lens. This phenomenon, 
known as the vignetting effect, 9 can be ignored if it is assumed that 
the lens is large enough so that the sum of the squares of the argu- 
ments of P L is always less than the square of the radius of the lens 
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for all values of x , x f , y , and y, for which U„ and T f are nonzero. 1 
With this assumption P, 4 may be taken as one. 

To further simplify equation (6) we introduce the notation 



5*. 

Xi " M ' 



y< ■ = ff > x = Xo + Xi , y = y + pi 



Then, neglecting the vignetting effect, 

U,(.v, , V.) = &i JJ U (x - £ ,. , y - y,)F[T f ] dx dy 



where 



F[T r ] = //_" T,(x, , y,) exp 



j2irM 

\d 



jdx, 



(xx f + yy f ) dx, dy, . (7) 



We observe that the disturbance in the image plane can be written 
as a convolution of the disturbance in the object plane and the Fourier 
transform of the transmission function of the phase plate evaluated 
at the spatial frequencies Mx/Xd, My/Xd. 

Let us now examine the form of F[T f ] more carefully. The transmis- 
sion function of the phase plate, from (4), is given by 



with 



and 



T,(r,) = (ae ia - 1) circ (r,/R ,) + circ (r,/R fo ) 

- (x 2 , + »})» 

1 r, SR 



circ (r,/R) = 



r, > R. 



It follows then that 



F[T,] = 



where 



M 2 



(%rMR° f 

, ia i.iVR , 'V Ad 



./, 



f'lirMR, 



+ 



MR, 'V \d 



\d 



(8) 



r = (x- + y 2 )\ 
and J, is the Bessel function of the first kind of order one. 



t Recall U. = for x 2 + y, 2 > #. , and T t = for x f 2 + y,z > R fo . 
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Now, at optical wavelengths, R f0 /\ is exceedingly large, and for all 
practical purposes, the second term in (8) can be approximated by a 
two-dimensional delta function with little error. 1 On the other hand, R° 
is typically of the order of the radius of the Airy disk of the diffraction 
pattern of the entrance aperture so that the first term cannot be so 
simplified. To emphasize this fact, R° will be written as 



R" f = y 



Xf 
ft- 



where 7 is a constant of the order of unity. (7 « 0.61 if R° is the radius 
of the Airy disk.) 

With this notation and the definition of M, equation (8) becomes 



F[T f ] 



\ 2 d 2 
M 2 



j(^r 



(9) 



Returning to the expression for the complex disturbance in the image 
plane (7), we conclude that 



Ufa , Vi) 



&'-hJL v (*-m>*-hI 



M 



./,i ';*</ + v*i\ 



■■[ 



R„ 



dxdy+ C/.I-17, - 



x, 



Mi 

M ' M 



(10) 



' (x 2 + y 2 )» 

Equation (10) is the fundamental result that will be employed through- 
out the remainder of this paper. 

III. SMALL PHASE OBJECTS-THE PHASE CONTRAST MICROSCOPE 

The form of equation (10) indicates that, in general, a complicated 
relationship exists between the object and the image. We will now 
demonstrate, however, that the simple theory, equations (2) and (3), 
may be employed if the ratio of the area of the phase object to the 
area of the entrance aperture is sufficiently small. This is usually the 
case for a microscopic system. 

We will assume that the phase object is centered in the entrance 



t We shall observe later that this approximation yields sharp jumps in intensity 
when the phase object possesses phase discontinuities. In practice some smooth- 
ing of the discontinuities in the intensity pattern will always result. (See Ref. 
10 for example.) We are neglecting these lower order effects for simplicity. 
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aperture and write U (x„ , y ) as 

U.(x. , y.) = circ [ (xl + yl) * ~\ + (e 1 * - 1)P. , 

where P has a value of one over the region covered by the phase 
object and a value of zero elsewhere. In the limit of a veiy small 
phase object, 1 the contribution of the second term in U to the integral 
in (10) approaches zero, and the integral can be approximated by 



_ ± rr eir C r«* - * y + (,J ~ tPH Lfi - f J 



dxdy. 



The convolution is most easily evaluated by first taking the Fourier 
transform of I and then taking the inverse transform of the results. 
These two operations yield 

I - £'" Ut)Jo(^t)dt (11) 

where ?-,• is the radial coordinate in the image plane. 

If we now restrict the coordinates in the image plane to the region 
corresponding to the image of the phase object, the ratio r\/MR will 
be small, and I can be approximated by 



I 



/,(<) dt = 1 - J n (2iry). 



Finally, with an appropriate choice of the radius of the central spot 
of the phase plate (and thus y ) , the contribution from the Bessel 
function J (2-y) can be made to vanish, and we conclude that, for a 
sufficiently small phase object and a suitable y, 

Ufa , yd w fl [ (ac ' a " 1} + cxp D*(~!/ ' ~M 
The intensity of the image is then 



I(xi , yd - I U t | 2 w jp 



(ac' a — 1) + exp 



W-l ■ -S)] 



which agrees with equation (3). The simple theory of the phase con- 
trast microscope is, therefore, a limiting case of the more general 
approach. 



f Consider a single cell under a microscope, for example. 
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IV. MACROSCOPIC SYSTEMS 

In many potential applications of phase contrast imaging the mag- 
nification of the system may be of the order of unity, and the area of 
the phase object may be comparable to the area of the entrance 
aperture. In these cases, the approximations that were employed in the 
last section cannot be used, and the complex disturbance in the image 
plane can, in general, only be determined via a direct application 
of equation (10). In order to provide some representative results, we 
have chosen a few simple, yet important, examples for which the 
integral in (10) can be easily evaluated. 

4.1 Background 

The purpose of any phase contrast system is to give a visual repre- 
sentation of the phase distribution in the entrance aperture. If a system 
is to yield meaningful results, therefore, the intensity of the image of 
the aperture should be essentially constant when no phase objects are 
present. In this case 

I7.fr. , V.) = eirc [&%&] 
and from (10) and (11) it follows that 



'« - w 



^'-i)[Mt)Jo(^t)dt + i 



(r, ^ MR.) 

(12) 

where the notation z = 2iry has been introduced. 

For our later purposes, it is convenient to define a normalized in- 
tensity, 7 = M 2 I. Plots of this function are shown in Figs. 2-5 for 
various values of z. Clearly, the response of the system to the illumi- 
nating light is critically dependent upon the size of the central spot of 
the phase plate. Where a = 1 and a = ir/2, a typical case that is em- 
ployed in phase contrast imaging, near uniform backgrounds result 
for values of z corresponding to spot sizes on the order of one-half the 
Airy disk or smaller.* For larger z, the background varies considerably 
until several rings in the diffraction pattern of the entrance aperture are 
covered by the phase spot.* We shall see later that these variations 
produce distortions in the intensity pattern of phase objects when they 



t The task of making a phase plate with a spot size this small may be difficult. 

* The phase contrast system will begin to behave like a direct imaging system 

for spot sizes larger than those shown, and thus large z values are not practical. 
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Fig. 2 — Background intensity distributions for spot sizes less than or equal to 
the Airy disk (« =1, a = jt/2). 



are placed in the aperture. The case a = 1, a = t/2, z = 1.5 will be of 
special interest to us, since in this case the integral in (12) is approxi- 
mately 1/2 except when r, is near the edge of the image of the aperture, 
the background is essentially constant, and several interesting results 
can be obtained. 

Frequently, a phase contrast system will employ a partially absorbing 
phase plate (a < 1) to improve the contrast between the image of the 
phase object and the background [see equation (3)]. Background in- 
tensities for a = 0.1, a = ir/2 and various values of z are shown in 
Figs. 4 and 5. There it may be observed that, in general, the back- 
ground intensity is reduced, although a bright ring appears in the 
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Fig. 3 — Background intensity distributions for spot sizes larger than the Airy 
disk (a = 1, a = ir/2). 



region close to the edge of the image of the aperture. (This ring has been 
observed experimentally by Kolettis.) If this edge effect is ignored, near 
uniform backgrounds are obtained for values of z that include five or 
six Airy rings. 

Taken together, the above results indicate that care must be exercised 
in designing a phase contrast system if near uniform backgrounds are 
to be achieved. 

4.2 Phase Disk 

Consider the "phase disk" defined by 

^.W.) = (*■«- 1) circ [^] + <4^] 

(0 ^ Rl l) ^ R.) (13) 

where </> is a constant. The geometry corresponds to a disk with radius 
R ( l) centered in the entrance aperture which introduces a constant 
phase shift 4> in the incident wave. Substituting in equation (10) and 
using the technique that was employed to derive equation (11), we 
obtain the intensity distribution 



?(r<) = ■ 
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(ae ia - IW* - l)I (l) (r t ) + I(r,)] + e<* | 2 , 

| (ae ia - 1)[(«'* - l)I (1) (r<) + I(r,)] + 1 | 2 , 

il/fl* 1 ' < r, g M# . 
In this expression I(r.) is the integral given in (11) and 

>(1) n2ry 



I U) (r,) = 



Ra 



£" y /.(^V ( : . i) >ii. 



MR„ 



1953 



(14) 



(15) 



Note that the integral I only takes into account the finite extent of the 

entrance aperture, while I (1) includes the size of the phase object as well. 

Plots of (14) are shown in Figs. 6-8 f or a = 1, a = ir/2 and z = 1.5. 

The value of z has been chosen to correspond to a near uniform back- 
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Fig. 4 — Background intensity distributions for spot sizes less than or equal to 
the Airy disk (a = 0.1, a = n/2). 



1954 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 

0.30 



0.25- 

0.20 

<M 

^ 0.15 
II 

0.10 
0.05 




0.4 0.6 

Pl/MRo 



Fig. 5— Background intensity distributions for spot sizes larger than the Airy 
disk (a = 0.1, a = w/2). 
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Fig. 6— Response to a phase disk (<£ = 0.1, R„W/R° = 0.1). 
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Fig. 7— Response to a phase disk (<t> = 0.1, R.M/R. = 0.5). 

ground when R ( „ x) = i2„ or when R"' = 0. It is clear from the figures 
that, for the small value of 4> chosen, the qualitative behavior of the 
phase distribution in the entrance aperture is reproduced quite nicely. 
The phase disk can be readily observed and the magnitude of the dis- 
continuity in I is on the order of #. t 

In Fig. 9 we have indicated the response of the same system with 
the exception that # has now been assigned a value of w. Note that in 
this case, the jump in 7 depends critically upon the size of the phase 
object. Moreover, it appears that a value of R { n exists for which no 
contrast results and beyond which the contrast is reversed. This sug- 
gests that for larger </> values, the system performs quite poorly. 

To explain this behavior, we need only to return to equation (14) 
from which it immediately follows that the magnitude of the dis- 
continuity in I is 

A/ = 2 Re {(ae ia - l)(e'* - l)[(e" - l)I a) (J*Bi n ) + I(MR™)]) 

= 2(a cos a - 1)(1 - cos#)(2I (l) - I) + 2al sin asin tf>. (16) 

For small <j>, 

A7 ;^' (2al sin a)<i>, 

while for larger <£, the full equation (16) must be employed. Recalling, 
that for z = 1.5, I .^ 0.5 for most values of r, , we conclude that, for 



t The increase in intensity in the outer portion of the image is due to the slight 
nonuniformity of the background in this region. 



1956 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 



;w 



0.60 












a = 


.0 






a = 


77/2 


0.55 


- 


z = 


.5 


0.50 








0.45 








0.40 


1 1 1 1 







0.2 0.4 0.6 0.8 

n/MRo 



Fig. 8— Response to a phase disk (<f> = 0.1, R.W/R. = 0.8). 
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Fig. 9 — Response to phase disks of various sizes (<£ = tt). 
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a = 1, a — if/2, z = 1.5, and small, 

as observed. 

When <f> = ir, however, equation (16) reduces to 

A7 = 4(a cos a - l)[2I lli (MR™) - I(MR™)]. 

From the definitions of the integrals (11) and (15), it is clear that for 
= 7T there exists a value of R l l) for which Al = and beyond which 
the contrast is reversed. 

In the preceding section, we observed that nonuniform backgrounds 
can result when the size of the central spot of the phase plate is not 
carefully chosen. Figure 10 demonstrates the response of a phase con- 
trast system to a phase disk when this is the case. The value of z has 
been chosen so that the spot size corresponds to the Airy disk of the 
diffraction pattern of the entrance aperture. The background corre- 
sponding to this value of z is shown in Fig. 2. Note that the phase 
disk simply perturbs the background pattern yielding an intensity pro- 
file that does not have the desired step function shape. 

4.3 Phase Rings 

The analytical results presented in Section 4.2 indicate that a dis- 
continuity in the radial coordinate of the phase distribution of the en- 
trance aperture can be qualitatively reproduced in the intensity pattern 
providing a sufficiently uniform background is chosen and the magni- 
tude of the discontinuity is small. In other cases, such as a large dis- 
continuity or a nonuniform background, rather poor results are obtained. 
We will now generalize these results to a more complicated radial phase 
distribution. 

Suppose the phase distribution in the entrance aperture consists of a 
central disk and a sequence of annular rings in each of which the phase 
is some constant 4> k . Let the radii defining the location of the phase 
discontinuities beR l „ k \ h =» 1, • • ■ , n — 1, and Ieti2J n) = R„ . Then U„ 
can be written as 

U (x a , iO - t, Ul k) circ 



Ri k> 



q 



(17) 



where 



uW = jexp [#>*.] - exp [#> A . Tl ] k = 1, • • • , n - 1 

[exp \j(j>„\ k = n. 
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Fig. 10— Response to a phase disk when the phase plate spot size equals the Airy 
disk of the entrance aperture (<£ = 0.1, R.W/Ro = 0.2). 

The normalized intensity distribution immediately follows as 
(ae ia - 1) £ Ui k) I ih \r t ) + exp [ft.] 

^ Ti £ MR™ 
(ae ia - 1) £ U™I m (rd + exp [ft l+1 ] 
Jlffli" <r, ^ MBS ,+1) , 1= 1, ■■■ ,n- 1. 
In this expression I (M (r«) is the generalization of (15) which is given by 

im w = f />(? Mfe')* (19) 

Plots of (18) are shown in Figs. 11-13 for a = 1, a = t/2, and z = 1.5. 
The assumed phase distributions in the entrance aperture are given at 
the top of the figures and the resulting intensity profiles are shown 
below. A qualitative reproduction of the phase distributions results in 

each case. 

The behavior of (18) in the limit of small phases, <f> k , explains some of 
the success that was achieved in the figures. From the definition of the 



PHASE CONTRAST IMAGING 



1959 



coefficients U[ k) it follows that, if all the <t> k are small and, in particular, 
if 0„ = 0, 



U? 



\K4>k - 0* +0 k = 1, • •• , n - 1 
I 1 k = n. 



Therefore, the jumps in intensity at the points r, = MR™, to first order 
in the phase jumps, is just 

Al(MRl k) ) fsa 2 Re {(«" - l)j(<fc +1 - 0*)I (n, (M/Ol 

= 2a sin «(</>, - 0, + l )I < "'(M/el*" ) ). (20) 

Now, in the cases shown in the figures, a = 1, a = ir/2, and 2 = 1.5. 
Recalling that for z = 1.5, 1 M « 0.5 we conclude that 

Thus, if the <£ fc are sufficiently small, and an appropriate value of z 
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Fig. 11 — Response to phase rings with decreasing phase. 
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Fig. 12— Response to phase rings with increasing phase. 

is chosen, the jumps in intensity are of the order of the jumps in phase. 

Figures 11-13 also indicate that, for the values of the parameters 
chosen, the piecewise constant behavior of the phase distributions is 
essentially reproduced in the intensity profile. The success here is 
mainly due to the fact that not only was z chosen to yield a nearly 
uniform background, but also the z value selected was relatively small. 
Under these conditions the dependence of the integrals I (W (n) is 
suppressed,* and the effect of the phase object is to introduce essen- 
tially piecewise constant perturbations to a nearly uniform background. 

To demonstrate what effects can occur for a larger z value, we have 
plotted in Fig. 14 the response of a phase contrast system to a 
sequence of phase rings when a = 0.1, a = ir/2, and z = 21.21 163. * 
This is a practical case to consider since small values of a are fre- 



t Consider equation (19) in the limit of small z. , 

t This value of z corresponds to a phase plate spot size that includes about 
five and one-half Airy rings. 
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qucntly employed in phase constant microscopy to improve the con- 
trast between the phase object and the background. (For our purposes, 
large z values are needed to give a fairly uniform background over a 
large portion of the image of the entrance aperture when a is small.) 
Unfortunately, the results that are observed are then quite poor. The 
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Fig. 13 — Response to phase rings with alternating phase. 

problem here is that the integrals I"' 1 (r,) are now highly nonlinear 
functions of r t , and, therefore, nonlinear perturbations of the order of 
the background are introduced. 1 

4.4 Semi-circular Phase Disks 

The results of the last section demonstrate that with an appro- 
priate choice of the parameters a, a, and z, circularly symmetric phase 



• Figure 5 demonstrates that for the values of the parameters chosen, the back- 
ground is of the order of a 1 = 0.01. However, from equation (20), we observe that 
the jumps in 7, which we may take as indications of the magnitudes of the perturba- 
tions of the background, are of the same order when <£t — <£t + i is small. 
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Fig. 14 — Response to phase rings (a = 0.1, z = 21.21163). 



distributions possessing discontinuities in the radial coordinate can be 
qualitatively reproduced in the intensity patterns. As a final example, 
we now examine the response of a phase contrast system to a simple 
angular discontinuity. 

Let the coordinates in the object plane be r„ and 0„ and consider the 
complex disturbance 



U (r , e B ) = circ (r /R ) exp [je(6 )}, 



where 



e(0.) = 



7i 



o ^ e„ ^ t 



7a 7T 



< e a < % 



(21) 
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witli yi and y 2 assumed constant. With this phase object, the entrance 
aperture is divided into two semi-circular regions in each of which 
the phase is constant. 

In order to calculate the intensity pattern that is produced by (21), 
we must compute the convolution integral (see equation (10)): 



R r 



(22) 



As in our previous analysis, we will accomplish this task by first taking 
the Fourier transform of I s and then taking the inverse transform of 
the result. 

Now, it is a well-known result in the theory of Fourier transforms 
that if a function cjir, 6) is separable in r and 6, i.e., g(r, 9) = 
Qr(r)ge{6), then the Fourier transform of g can be expressed by the 
following infinite series of Hankel transforms: 



where 



and 



F[g] = £ c k (-]) k e ik *K k [9r(r)] 



dO 



II 



\g r (r)] - 2tt f rg r (r)J k (2irrp) dr. 

•In 



In these expressions p and <f> are the coordinates in the transform 
space, and J k is the Bessel function of the first kind of order k. 
Applying this result to (22) yields 



where 



F[I.] = circ (^ E ft<-j)V** f " rJ k (2*rp) dr (23) f 



c k = 



Kexp fr/i] + exp [jy 2 ]), k = 

k even . 

JS (exp [jv 2 ] ~ exp [jy t ]), k odd 



(24) 



t Here we have used the easily established result that 

F\ -£-/,^ / I = circ (pRo/y). 



' ^Ji^yr/R,,) 
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The appropriately normalized inverse transform of (23) is 

I„ = f f F[I.]p exp [j2irf iP cos (0. - <*>)] dp d<f> 
Jo Jo 

oo ny/Ro /•«» 

= 4tt 2 2 c, exp [jkdi] / p / rJuClirrp) drJ k (2nr iP ) dp. 

Here 0, is the angular coordinate in the image plane, and we have already 
introduced r\ = 1\/M for later simplicity. Now, letting t = 2irR p, 
s = r/R c , and substituting expressions (24) for the coefficients c k 
finally gives 

1 2 

I„ = -(exp [jy l ] + exp [jy 2 ])h (ri) + - (exp [jy 2 ] - exp [jy,]) 



todd 

where 



Z k*o ^ ( 25 ) 



fc fc ( r .) = J'' T | £ s / fc (rf) dsJ k {^) dt. (26) 

Note that the integral h (ri) is just the integral given in (11), and that 
the integrals A§(r<) are higher order generalizations of it. 

The intensity in the image plane is obtained from equations (10) 
and (25) as 

(| (ae ia - 1)1. + exp [jy 2 ] | 2 , £ 0< £ *■ 

/(r,- , 8<) = | ^ '< ^ Mtf . 

l| (ae' a - 1)1. + exp [jyA \ 2 , * < 6 { < 2w (27) 

Clearly 1 depends upon both the radial and the angular coordinates in 
the image plane while the phase distribution over the entrance aperture 
depends only upon 0. Therefore, one would not expect, in general, that 
the semicircles would be visible in the intensity pattern. We will now 
demonstrate, however, that if a small value of z = 2iry is used for 
which the background is essentially uniform, and if 72 — 7i is small, 
the semicircles can be reproduced, and the jump in the normalized 
intensity will be proportional to 72 - 7i , at least to first order. More- 
over, as in our previous analysis, if a = 1, a — ir/2, and z = 1.5, the 
first condition will be satisfied, and the proportionality factor will be 
approximately unity. 

To observe these facts we need only to substitute the first term of the 
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power series for the Bessel functions in (26) and conclude that for z 
sufficiently small 



kMRJ (Tk\y(k + 2)(2fc + 2) 

Thus for small z, the coefficients h k (r t ) are small, and they decrease 
extremely rapidly as k increases. Hence, if y 2 — 7, is sufficiently small, 
the first term in (25) will dominate the second and I, can be approxi- 
mated by 

I. « £(exp [JYi] + exp [frj)ft.(r<). 

With this approximation the normalized intensity is approximately 

|(ae' a - l)(exp [jy t ] + exp \jy»])h.(r t ) + exp [jy a ] | 2 , 

£ di £ir 

I(r< , 6i)tt\ ^ ?-, ^ Mfl. . 

Kae'° - l)(exp [jy,] + exp tfy 2 ])fc (r,) + exp [jy,] | 2 , 

f < J, < 2ff. 

Now, if the background is essentially uniform, the integral h e (r { ) in 
the above expression is essentially constant, and we observe that, for all 
practical purposes, I is constant in each of the semicircular regions 
5^ 0, ^ 7T, 7T < 0, < 2ir. Furthermore, the jump in 7 between these 
two regions is just 

Al « 2 Re {§(ae' a - l)(exp [-772] - exp [-jy.])(exp [jy,] 
+ exp [jya])h,(ri)} 

= 2a sin a /t„(r,) sin (7^ — 7,) 

~ 2a sin a ft (r,)(7 2 — 7,) . 

From our previous work, we may recall that, for z = 1.5, Mr t ) ~ 0.5, 
and thus, for a = 1, « = tt/2, and 2 = 1.5, 

A7 « (7 2 - 7,). 

The same conditions that permitted radial discontinuities in the phase 
distribution to be qualitatively observed in the intensity pattern 
appear to permit a simple angular discontinuity such as (21) also to 
be observed. The representative results of numerical computations of 
the exact expressions (25), (26), (27) shown in Fig. 15 bear this out. 
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Fig. 15 — Response to a semi-circular phase disk. 
V. CONCLUSIONS 

A general expression for the intensity distribution that is produced 
by a phase contrast imaging system with a circular phase plate was 
derived. It was shown that the results reduced to the well known 
expression for the phase contrast microscope when the size of the phase 
object is small compared to the area of the entrance aperture of the 
system. To obtain the intensity distribution for larger phase objects, 
a convolution integral must be evaluated by analytical and/or nu- 
merical techniques. 

The intensity patterns that are produced by phase disks, phase 
rings, and semicircular phase distributions were derived, and the re- 
sults of numerical computations that were based on these derivations 
were studied. In general, it appeared that qualitative reproductions 
of these simple phase distributions could be observed in the intensity 
patterns if the size of the phase spot on the phase plate was chosen 
to yield a uniform background, if the resulting parameter, z, was 
small, and if the magnitudes of the phase perturbations were small 
enough so that the approximation e>* ~ 1 + jc/> held. Moreover, it 
appeared that, if these conditions were not met, rather poor reproduc- 
tions of even these simple phase distributions could result. 
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